Chimera states on a flat torus 
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Discovered numerically by Kuramoto and Battogtokh in 2002, chimera states are spatiotemporal 
patterns in which regions of coherence and incoherence coexist. These mathematical oddities were 
recently reproduced in a laboratory setting sparking a flurry of interest in their properties. Here we 
use asymptotic methods to derive the conditions under which two-dimensional chimeras, similar to 
those observed in the experiments, can appear in a periodic space. We also use numerical integration 
to explore the dynamics of these chimeras and determine which are dynamically stable. 
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In nature, arrays of oscillators often synchronize and 
oscillate with a single frequency. This phenomenon can 
be observed in diverse systems ranging from laser ar- 
rays [1 and Josephson junctions [2H3], to populations 
of fireflies [5] and heart cells [6 . While incoherence and 
synchronization are ubiquitous in these systems, other 
complex patterns are also possible. 

The standard mathematical paradigm for modeling ar- 
rays of coupled oscillators is the Kuramoto model [7]. 
Kuramoto approximated the complex Ginzburg-Landau 
equation and demonstrated that under weak coupling, 
amplitude changes can be ignored |8|. Thus, oscillators 
can be approximated as coupled only through their phase 
with dynamics governed by 
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where fa is the phase of oscillator z, uji is the natural 
frequency, K is the coupling strength, and a is the cou- 
pling lag [9]. For arrays with narrow unimodal natural 
frequency distributions and no lag (a = 0), a first or- 
der phase transition occurs. Below a critical coupling 
strength, oscillators remain incoherent and above this 
threshold, they begin to synchronize [TQ| [TTj . 

In 2002, Kuramoto and Battogtokh observed that with 
nonzero lag and nonlocal coupling, surprisingly, regions 
of coherence and incoherence can coexist even for iden- 
tical oscillators (u)i = uu\/ i) [12] . Abrams and Strogatz 
described this hybrid state as a "chimera" [13]. Since 
then, chimera states have been observed in various sys- 
tems including: two groups of oscillators with no spatial 
extent p~4[ [15] , a ring of oscillators [12] [13] [16] , an infinite 
plane [T7] [18] , and a periodic two-dimensional space [19] . 
These states display a variety of spatial patterns includ- 
ing stripes, spots (see Fig. [I]), and spirals. 

Recently, two experiments observed chimeras in lab- 
oratory settings for the first time [20 . Tinsley, Nkomo 
and Showalter used photo-excitatory feedback to couple 
two populations of discrete chemical oscillators [21 . Oc- 
casionally, one population synchronized while the other 



remained incoherent. This is consistent with the analyt- 
ical results in ref. [T4] . 

Meanwhile, Hagerstrom et al. used a computer with 
feedback from a camera to control the phase modula- 
tion induced by a spatial light modulator [22]. This cre- 
ated a physical realization of a two-dimensional iterated 
map with nonlocal coupling and periodic boundary con- 
ditions. They observed the formation of chimera states 
as one-dimensional stripes of incoherence. Omel'chenko 
et al. studied a similar system and produced both sym- 
metric spot and stripe patterns in numerical experi- 
ments [19] . These patterns have yet to be explained from 
an analytical perspective. 

Analysis — To determine conditions for the existence 
of these chimera states, we examine a two-dimensional 
array of oscillators in a space with periodic boundaries: 

T 2 = S 1 x S 1 = {{u,v)\u e [0, 27r), v £ [0,2tt)}. 

This can be interpreted as the surface of a torus where 
the u- and ^-coordinates correspond to the toroidal and 
poloidal angles (we disregard any effects of surface cur- 
vature) . 

We consider a generalization of the traditional Ku- 
ramoto model to a continuous distribution of oscillators: 

=w-J a G(r, r') sin(0(r, t) - 0(r', t) + a)dr' , 

(1) 

where G(r, r') is a continuous coupling kernel. Following 
the approach of Kuramoto and Battogtokh [12 , we shift 
into a rotating frame with angular frequency Q (to be 
determined later) and define a complex order parameter, 

i?(r,t)e i0(r '^ = / G(r,rV^ r '' £ W, (2) 

J T 2 

resulting in a new governing equation: 

= A - R(r, t) sin(0(r, t) - 9(r, t) + a), (3) 
where = (j) — Vtt and A = uo — Vt. 
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FIG. 1. Stripe and symmetric spot chimeras. Each panel contains a 200 by 200 grid of oscillators. Each pixel corresponds 
to a single oscillator with color representing the phase. Panels A and B display a stripe chimera state with a\ — 0.2322, 
Ci2 — 0, 5 = 0.1065, while C and D display a symmetric chimera state with a coherent spot and with a\ — 0.1366, a<z — 0.1366, 
5 = 0.1494. Numerical integration indicates that these states are stable. 



In regions where R(r) > |A|, stationary solutions to 
Eq. (|3| exist. Oscillators become phase- locked and ro- 
tate at a fixed angular frequency. In regions where 
R(r) < |A|, stationary solutions are not possible. In- 
stead, oscillators drift at a nonzero phase velocity and 
satisfy a stationary probability density. Kuramoto and 
Battoogtokh observed that Eq. ([3| can be used to elim- 
inate e z6> ( r ' £ ) from Eq. (J2|. This yields a self-consistency 
equation for solutions to Eq. (J3|: 
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This functional self-consistency equation is effectively 
infinite-dimensional and yields little insight. To proceed, 
we define a simple kernel representing nonlocal coupling: 

G(r, r') = /o \ [1 + k(cos(u - u') + cos(v - v'))} . 
( 27r ) 

This is the leading order approximation to the two- 
dimensional von Mises distribution, the circular analogue 
of the Gaussian, and can also be interpreted as a pertur- 
bation off of all-to-all coupling for k <C 1. This choice 
allows us to remove explicit dependence on u and v from 
the integrals and express the order parameter as follows: 

R(r) e i@ ( r ) = c + di cos(u) + d 2 cos(v), (5) 



where (/(r')) = j^y J T 2 f(*')dr' and 

c = e i/3 (/i(r / )e i0(r/) ) (6a) 
di = «e i/3 (h(rV e(r,) cos(i/)) (6b) 
d 2 = Ke ip (h(r f )e i@ ^ cos(v')). (6c) 

Without loss of generality, we define = at the 
point (uo,vo) = (§,§), resulting in real- valued c. We 
use Eq. ([5| to eliminate R and 6 from system (J6|, thus 
reducing infinite dimensional functional Eq. Q to a set 
of six algebraic equations (3 real and 3 complex) with six 
variables (c, Re(di), Im(di), Re^), Im(^2) and A): 



A 2 — \c + d\ cos(u') + c?2 cos(v')\ z 
+ d* cos(u f ) + cos(v') 



(7a) 



/ A - ^ - \c + dx cos( U ') + d 2 co S (v')\ 2 > 



c* + d{ cos(u') + dZ, cos(u') 



(7b) 



/ A - ^ - \ c + d x cos(u') + d 2 cos(v')\ 2 \ 
di = Ke ' { — : — — cos(v ) ), 



c* + d\ cos(u') + c?2 cos(u') 



(7c) 
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where * denotes complex conjugation. 

System Q can be solved asymptotically for k = e <C 1. 
Motivated by results from one dimension [T6] . we make 
an ansatz for the scalings 



C ~ 1 + C\€ + C26 2 

di ~ (ai + i6i)e 2 
d 2 ~ (a 2 + ^2)e 2 
A - 1 + Aie + A 2 e 2 , 



substitute these expressions into Eq. (7a), and retain 
terms to 0(y/e): 



1 + O(e) = 1 + y/2y/A 1 -c lv ^ + ^(e). 

For this to be satisfied, Ai = c\ is required. Expanding 
system ([7]) to leading order we obtain 



ci 



i/3i — V2 y/ 5 — CL\ cos(u') — CL2 cos(V) ^ 



(9a) 



ai 



• Z&i = — COs(l/) y/S — CLi COs(u f ) — CL2 COS(V) 



(9b) 



a 2 + ibi = —V2(^ cos(v f ) y/S — a\ cos(u') — CL2 cos(V ) 

(9c) 

where we have defined 5 = A2 — C2 for convenience. 

These three complex equations contain six variables: 
ci, ai, 61, a2, ^2 and (5, and one parameter: /?]_. Practi- 
cally, it is easiest to parameterize the system by 5, solve 




FIG. 2. Solutions to system [9] The inset displays the full set 
of solutions. The main panel shows the quadrant containing 
only ai , ci2 > 0. Loops correspond to symmetric spot chimeras 
(ci2 — ai) and two types of equivalent stripe chimeras (ai = 
and a2 = 0). Branches (a± / a2, both nonzero) correspond to 
two types of asymmetric spot chimeras, one with 5 > and 
one with S < 0. 
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FIG. 3. Cross sections of Fig. [2] The magenta (dash dot) 
curves represent unstable chimeras, the blue (solid) represent 
stable chimeras, the green (dashed) correspond to modulated 
drift states, and the red (solid, a\ — 0) correspond to uniform 
drift states. Panel A displays the stripe chimera loop {ai — 0). 
Panel B displays the symmetric spot loop (a\ — a^). 



for a\ and CL2 using the real parts of implicit Eqs. (9b) 
and (9c), and then deduce the remaining values. This 



yields a set of solutions corresponding to various types of 
chimera and drifting states. 

Results — Fig. [2] shows solutions to system It con- 
tains three closed loops connected by two branches. The 
loops lie in the planes a\ = 0, d2 = 0, and a\ = CL2 
whereas the branches have a\ 7^ CL2 (both nonzero). 
Because Eq. ([5| is invariant under the transformations 
(di,u) — » {—d\,u — 7r), (d2,v) — > (—d2,v — 7r), and 
(di : u) —> (d,2,v), the quadrant with nonnegative a\ and 
CL2 contains all of the distinct solutions to system ([9|, and 
the loops satisfying a\ = and CL2 = are essentially the 
same. 

The chimera states these solutions describe have locked 
and drifting regions separated by a contour where R = 
|A|. Expressing this boundary in terms of the variables 
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FIG. 4. Asymmetric spot chimeras. Panels A and B display an asymmetric chimera from the lower branch with a\ — 0.07635, 
d2 = 0.02929, 5 = —0.02608, while C and D display an asymmetric chimera state from the upper branch with a\ — 0.2089, 
d2 = 0.04844, S = 0.1769. Numerical integration reveals that these states are unstable. 



from Eq. ^ yields 



S = a\ cos('u) + d2 cos(v). 



(10) 



This boundary takes on various forms depending on the 
values of a±, ct2 and S. 

When ai = or 02 = 0, system ^ reduces to the 
equations analysed in ref. [16 j. The solutions along this 
branch represent one-dimensional stripe chimeras — the 
order parameter varies in only one spatial dimension. 
One such chimera is displayed in Figs. [ljA) and[ljB). 
Fig. [3jA) describes the various solutions that are found 
along this branch and indicates their stability. These 
chimeras are qualitatively similar to those observed in 
experiments [22] . 

When a\ = <22, the boundary is symmetric about the 
line u — v and resembles a circle or a square. These 
chimeras consist of an incoherent spot surrounded by a 
coherent region or a coherent spot surrounded by an in- 
coherent region, similar to those studied in ref. [19 . An 
example of such a state can be found in Figs. [lJC) and 
[TjD). Solutions along this loop are described in Fig.j3^B). 

When a\ ^ <22, the boundary resembles an ellipse or 
a rhombus, both of which are asymmetric (they are not 
invariant under reflections about the line u — v). There 
are two such branches, one with 5 > and one with 
S < 0. These solutions also have regions of incoher- 
ence surrounding or surrounded by coherence. Examples 



of these previously undiscovered chimeras are found in 
Fig.g 

We determined the stability of predicted chimera states 
by numerically integrating Eq. Qwith initial conditions 
determined by our theory. Fig7[5] plots the fraction of 
oscillators in the drifting region and the asymmetry of 
that region (defined as min(|ai/a2| , ^2/^1 1)) as a func- 
tion of Pi. The stripe chimera loop contains both stable 
and unstable domains, and is nearly identical to the so- 
lution curve described in ref. [16 (see Fig. ^A)). The 
symmetric spot loop also contains both stable and un- 
stable domains (see Fig. [3^B)) and is similar in shape 
to the stripe chimera loop, but it has differing stabil- 
ity regions. Both asymmetric branches are unstable, and 
nearby states evolve within planes of fixed j3\ to solutions 
along the stripe or symmetric spot loops. 

Conclusions — This work reveals a new type of chimera 
state in a two-dimensional periodic space. These asym- 
metric chimeras are unstable and appear in simulations 
only as transients. Nonetheless, understanding where 
these asymmetric chimeras appear is essential to under- 
standing where stripe and spot chimeras are stable. Sta- 
ble chimeras are created via continuous bifurcation off of 
modulated drift states, while unstable chimeras appear 
via continuous bifurcation off of the fully synchronized 
state. As j3\ increases, symmetric and stripe loops inter- 
sect with asymmetric branches resulting in unexpected 
changes in stability. Eventually, a stable chimera collides 
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equilibrium states. On a torus, single spirals are excluded 
by the periodic boundaries — they are topologically im- 
possible. However, on an infinite plane, finite-sized non- 
spiral chimeras have not yet been observed and may not 
exist. Networks of oscillators are often of interest, and 
do not necessarily reproduce the topological properties of 
any simple metric space. We hypothesize that, on arbi- 
trary networks of oscillators with more complex structure 
[23] , previously unobserved chimeras are possible. Inves- 
tigating these chimeras may shed light on breakdowns of 
synchrony observed in nature. 
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FIG. 5. Dynamics of the fraction drifting and asymme- 
try. Solid lines indicate stable states while dashed lines in- 
dicate unstable states. Red curve: symmetric spot chimeras, 
blue curve: stripe chimeras, green curve: asymmetric spot 
chimeras (5 > 0), purple curve: asymmetric spot chimeras 
(8 < 0). Unstable states evolve along planes of fixed f3\ to 
nearby stable states (with a different fraction drifting and/or 
asymmetry, see text for definition) as indicated by the arrows 
in panel A. Panel B contains a two-dimensional projection of 
panel A showing the fraction drifting as a function of fii for 
the various chimeras. 



with an unstable chimera causing both to be destroyed 
in a saddle-node bifurcation. 

Our analysis reveals the complex regions of parameter 
space in which two-dimensional chimeras reside. The in- 
tricacy of bifurcation diagrams [3] and [5] elucidate why it 
is difficult to reproduce chimera states in numerical sim- 
ulations and experiments. Stable chimeras only exist for 
narrow ranges of the system parameters. Consequently, 
small changes in the lag parameter a (or equivalently j3) , 
can cause stable chimeras to vanish or new chimeras to 
appear. 

These findings also highlight the impact of topology on 
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